function diff = jtj_num( delta, q)

    % Calculates job-to-job rate given grid for delta(m) and searcher distribution 

    G=(q -q(1))./(q(end)-q(1));
    g_p=G(2:end)-G(1:end-1);
    diff=g_p'*delta(2:end);

end
    